Optimal truss design with MOHO: A multi-objective optimization perspective

This research article presents the Multi-Objective Hippopotamus Optimizer (MOHO), a unique approach that excels in tackling complex structural optimization problems. The Hippopotamus Optimizer (HO) is a novel approach in meta-heuristic methodology that draws inspiration from the natural behaviour of hippos. The HO is built upon a trinary-phase model that incorporates mathematical representations of crucial aspects of Hippo’s behaviour, including their movements in aquatic environments, defense mechanisms against predators, and avoidance strategies. This conceptual framework forms the basis for developing the multi-objective (MO) variant MOHO, which was applied to optimize five well-known truss structures. Balancing safety precautions and size constraints concerning stresses on individual sections and constituent parts, these problems also involved competing objectives, such as reducing the weight of the structure and the maximum nodal displacement. The findings of six popular optimization methods were used to compare the results. Four industry-standard performance measures were used for this comparison and qualitative examination of the finest Pareto-front plots generated by each algorithm. The average values obtained by the Friedman rank test and comparison analysis unequivocally showed that MOHO outperformed other methods in resolving significant structure optimization problems quickly. In addition to finding and preserving more Pareto-optimal sets, the recommended algorithm produced excellent convergence and variance in the objective and decision fields. MOHO demonstrated its potential for navigating competing objectives through diversity analysis. Additionally, the swarm plots effectively visualize MOHO’s solution distribution of MOHO across iterations, highlighting its superior convergence behaviour. Consequently, MOHO exhibits promise as a valuable method for tackling complex multi-objective structure optimization issues.


Introduction
Engineers dealing with structural design often have competing objectives that can adversely affect each other.For instance, they might attempt to reduce the structure's weight for economic purposes, even if it would increase safety by strengthening the structure.MO optimization problems are characterized by multiple objectives, requiring mathematical optimization strategies to deal with them effectively.Multi-objective optimization problems typically need more than simple solutions.Instead, their responses usually manifest as a series of optimal solutions embodying a balance between the conflicting objectives being optimized.This tradeoff guarantees that achieving a particular purpose will inevitably result in the degradation of other objectives [1,2].The Pareto front provides the designer with various feasible alternatives by meeting all the constraints; among these, the final design can be selected.Applying metaheuristic optimization methods has become an increasingly popular research direction for multi-objective structural design.When dealing with difficult circumstances, metaheuristic algorithms have become highly sophisticated techniques concentrating on minimizing or maximizing an objective function to arrive at an optimal solution.These algorithms efficiently traverse the solution space by utilizing upper-level searching techniques.
Numerous researchers have introduced innovative optimization algorithms across various applications, prioritizing statistical validation and rigorous experimentation to enhance the convergence time and solution quality compared to existing methods.These algorithms address complex optimization problems in domains such as engineering design.[27], metaheuristic approaches, and medical diagnostics [28], demonstrating encouraging developments in effectively resolving practical issues [29] and different strategies [30,31], such as boundary update techniques, multi-hybrid approaches [32], multi-strategy infused metaheuristics [33][34][35].Each algorithm incorporates new metaheuristic methodologies to handle specific optimization goals and enhance overall performance.
These approximate optimization methods provide sufficiently good solutions for complex situations within a reasonable time frame.Population-based metaheuristics are particularly noteworthy since they are effective tools that may be used to solve MO optimization problems [36].Some notable multi-objective algorithms include the MO version of Atomic Orbital Search.[37], MO version of the Material Generation Algorithm [38], MO Crystal Structure Algorithm [39], MO Chaos Game Optimization [40], NSGA-2 [41], MO Arithmetic Optimization Algorithm [42], MO Thermal Exchange Optimization [36], MO Passing Vehicle Search Algorithm [43], MO symbiotic organism search [44], and MO heat transfer search [45].Many researchers have Improved MO optimization algorithms with unique approaches, such as decomposition-based MO heat transfer search.[46], improved MO particle swarm optimization [47], an indicator-based multi-SSM algorithm [48], MO improved marine predators algorithm [49], MO structural optimization using improved heat transfer search [50] Enhanced MO GWO with levy flight and mutation operators for feature selection [51], and a two-archive MO multi-verse optimizer for truss design [52].
The primary approach to estimating a Pareto front combines a non-dominated (ND) sorting technique with the population-based concept of meta-heuristics in MO optimization.The Pareto archive is continually refined by updating it with data from the current population and the archive of the previous iteration through the iterative reproduction of design solutions.This process was repeated until the Pareto front of the design problem was reached.Critical performance factors are optimizing the parameter settings using a self-adaptive method and preserving the search diversity through efficient clustering approaches.The main intention is to significantly increase the search intensity and variety because MO metaheuristics require extensive design space exploration while maintaining a high convergence rate, which is a more difficult challenge than single-objective metaheuristics.
The "No Free Lunch" (NFL) [53] theorem serves as a reminder that no individual metaheuristic is universally capable of solving all real-world problems.This understanding has catalyzed the advancement and enhancement of diverse metaheuristic methods.This research introduces an MO version of the Hippopotamus Optimization algorithm (HO) [54], drawing inspiration from the observed behaviors in hippos.Emphasizing traits and actions, such as their semi-aquatic lifestyle, herbivorous diet, and defensive strategies against predators, the algorithm incorporates their formidable jaws and warning vocalizations.Furthermore, it explores the adaptability of protective behaviors and social dynamics, aiming to understand their applicability in various contexts through insights from the hippopotamus behavior.
It is interesting to evaluate the performance of a newly created multi-objective metaheuristic in various engineering design challenges.This study applies MOHO to different five-truss structures: 10-bar truss, 25-bar truss, 60-bar ring truss, 72-bar truss, and 942-bar tower truss.The design challenge is decreasing the structure's mass while decreasing the maximum nodal displacement, a structural stiffness indicator.The above objectives were pursued within the limitations of area and stress.HO's mathematical model for organizing hippopotamuses within herds considers factors such as dominance and proximity to optimize their positioning in aquatic environments.defensive strategies for hippos against predators, including vocalizations and tactical movements, to ensure herd safety.Additionally, it outlines a behavior in which hippos escape predators by seeking refuge in water bodies, enhancing their survival prospects.Overall, this study explored methods to improve the safety and survival of hippopotamus herds through mathematical modeling and behavioral analysis.The primary progress of this investigation and its evolution, which surpasses the present contemporary, are outlined as follows.

A multi-objective version of the unique hippopotamus optimization algorithm (MOHO)
was applied to five planar and spatial truss structures to minimize the maximum nodal displacement and structural mass subjected to area and stress constraints.
3. Four commonly used performance metrics were used to assess the algorithms' efficacy statistically: the Hypervolume Index (HV), Generational Distance (GD), Inverted Generational Difference (IGD), and spacing-to-extent (STE).In addition, each algorithm's best Pareto-front plots were closely examined to evaluate the qualitative behavior.In addition, the algorithms were ranked for a thorough study using Friedman's test at the 95% significance level.
4. The results presented a fresh outlook on the benefits and drawbacks of MO optimization techniques in addressing diverse contradicting objectives in structural optimization.
The outline of the paper is as follows: ▪ The 2 nd section presents a description and mathematical modeling of the elementary Hippopotamus Optimization Algorithm (HO).
▪ 3 rd section elaborates on the proposed MOHO and outlines the formulations of MO structure optimization problems.
▪ In the 4 th section, an experimental assessment of the MOHO optimizer and a comparison with other prominent algorithms for addressing truss bar problems are presented.
▪ Section 5 discusses the performance metrics and enhanced Pareto fronts using diversity curves and swarm plots.
▪ Finally, Section 6 concludes the study and discusses future work to explore the capability of MO algorithms.

Hippopotamus optimization algorithm
The hippopotamus, a captivating African vertebrate mammal, thrives in semi-aquatic habitats, such as rivers and ponds, displaying social behaviors within pods.Despite the challenges in gender determination, their herbivorous diet and exploratory nature drive them to investigate alternative food sources [61].They are deemed one of the most dangerous mammals with immense strength and territorial behaviour, yet predators avoid confronting adult hippos because of their size.Defensive strategies include aggressive postures, loud vocalizations, and the rapid retreatment of water bodies.Inspired by the observed behavioral patterns, hippopotamus groups consist of various members, with calves prone to wandering and becoming prey [62].Defensive behaviors include rotation toward predators, employing powerful jaws, and fleeing toward water sources for safety [63].This algorithm uses two exploration stages and one exploitation phase after initial random solutions are produced.These phases are perceived to be superior to other algorithms for truss design problems.

Exploration phase-1 (position update in river or pond)
A herd of hippos consists of adult females, calves, adult males, and a dominant male leader.A value iteration process establishes dominance.While females encircle the dominant male, they defend the herd and its area.When a male reaches adulthood, the dominant male drives him out, forcing him to fight for supremacy elsewhere.Table 1 shows the mathematical positions of the hippos and their updates in the herd's habitat.Stage 1 in Table 1 represents the population of hippos after initialization within the upper and lower bounds of the lake or pond.Stage 2 shows five different scenarios based on random numbers h 1 and h 2 .Stages 3 and 4 demonstrate the position of immature hippos from their mothers within the herd, otherwise separated based on the T value.Finally, in Stage 5, the position is updated based on its objective value.

Exploration phase-2 (defense against predators)
In this phase, hippopotamuses find safety in their herds, deterring predators due to their large size and collective presence.However, young and sick individuals are more susceptible to predator attacks.When threatened, hippopotamuses emit loud vocalizations and may approach predators to avoid potential threats.Table 2 demonstrates a mathematical representation of half of the population of hippos with exploration phase 2.
In the second phase of exploration, which spans stages 6 to 11, the algorithm simulates the defense mechanisms of hippos against predators.
In Stage 6, the positions of the predators were determined within the search space.At stage 7, the algorithm evaluates the position of each predator relative to each hippo, triggering the corresponding reactions from the hippos.Stage 9 involves random movements of hippos, influenced by a Levy distribution, allowing search space exploration.Stages 10 and 11 determine whether predators hunt the hippos or if they successfully escape from predators.Overall, this phase simulates a dynamic interaction between predators and hippos, guiding the exploration process and helping to prevent the algorithm from becoming trapped in local minima.

Exploitation phase (escaping from the predator)
In this exploitation phase, as shown in Table 3, hippos typically run to the closest body of water for safety when separated from the herd and attacked by lions or spotted hyenas.This strategy improves local search capabilities by simulating the habit of seeking refuge nearby.This behavior involves creating a random site close to the hippopotamus's location to increase is the position of a male hippo.D do min ant-hippo is the position of the dominant hippo.i = candidate solutions.j = decision variables.
: Mgi = mean of some randomly selected hippos.
T represents the distance of an immature hippo from its mother.the cost function value.Iteratively, the hippopotamus adjusts its position to guarantee its proximity to safety.According to Table 3, during the exploitation phase, hippos can use the value of S from Stage 14 to determine a safe location near where they are right now to protect themselves from predators.Finally, a position update for the hippos was made using their objective values.

MO optimization definitions
In optimizing truss structures with multiple objectives, the focus is on simultaneously achieving optimization across various goals.While the term "multi-objective" typically addresses problems involving up to three objectives, the emergence of "many-objective optimization" many-objective optimization tackles the challenges posed by numerous objectives.Resolving conflicts between objectives is a significant challenge in multi-objective optimization, requiring specific approaches.The traditional relational comparisons between solutions are inadequate when there are several criteria.Therefore, alternative operators, such as the Pareto Defensive behavior of hippos based on the factor F predator j for protection against predators.
RL * is a random vector based on levy distribution.
. . .; N and j ¼ 1; 2; . . .; m If hippos have been hunted, other hippos will replace them in the herd; otherwise, the hunter will escape, and the hippos will return to the herd. https://doi.org/10.1371/journal.pone.0308474.t002 is the position of the hippo to find a closet safe place.
S leads to more suitable local search (exploitation).
r 12 is normally distributed random number. 15 https://doi.org/10.1371/journal.pone.0308474.t003 dominance operator, are essential to assess the relative superiority.From a mathematical point of view, Pareto optimality defines a set of solutions deemed nondominated and optimal for a given problem, particularly for MO optimization.These solutions constitute the MO optimization solution set, which thoroughly depicts realistic trade-offs between several objectives.This set, which illustrates the optimal solutions possible within the objective domain, is frequently represented and visualized as a Pareto optimal front (PF).Further clarification regarding the principles of domination and associated ideas can be obtained from Fig 1, which graphically presents these concepts.

Multi-Objective Hippopotamus Optimization algorithm (MOHO)
The MOHO algorithm starts with a random population (hippos) size N.It generates solutions X 1 , X 2 , . .., X N within the specified search space defined by the lower and upper bounds of the problem for each dimension (D).The fitness function evaluates and assigns a fitness score to each candidate solution, guiding the selection process toward the optimal solutions.The fitness of each candidate solution was assessed using the fitness function provided for all initial solutions.The main loop of the HO algorithm starts with the first iteration and continues up to a specified number of generations.The best candidate solution (Xbest) and its corresponding fitness (fbest) were updated for each iteration.The loop was divided into three stages (two exploration phases and one exploitation phase).Half of the total population explores the search space in exploration phase 1; each individual calculates two potential positions (X_P 1 and X_P 2 ) based on its current position, the best solution found so far (X best ), and a randomly selected mean group (MeanGroup).The selection of potential positions involves randomness, which is affected by various parameters.Individuals update their positions based on their potential positions if their fitness is improved.The remaining individuals in the population defend themselves against potential threats from predators.Each individual calculates a potential position (X_P 3 ) based on its current position, randomly generated predator position, and Levy flight.The decision to update positions is based on whether the new position improves fitness compared to the current position.All individuals perform a localized search to escape potential predators and exploit promising regions.Each individual calculates a potential position (X_P 4 ) based on its current position and randomly selected direction (D).Positions are updated if the new position improves fitness.The best solution found thus far, along with its corresponding fitness value, was stored for M objective functions as Pareto fronts.For every iteration, all fronts' non-domination (ND) sorting process helps achieve superior ND fronts.
Fig 2 shows a flow chart of MOHO with two exploration phases, dividing the total population by half size for phase 1 (position update in rivers or ponds) and half size for phase 2 (defense against predators).After passing through any of the exploration phases, phase 3 (escaping from predators), the exploitation phase, will be activated, and the positions of hippos will be updated accordingly.The evaluated solution at the end of this phase is the non-dominated Pareto optimal front.

Formulation of the truss design problem
MO truss design problems aim to determine the best design variables to minimize the structure's weight and optimize the nodal displacement while adhering to the tensile and compressive stress constraints.The computational problem is defined as follows.
to minimize the structure mass Subject to: Here, Ax is a vector for the cross-sectional area; ρ i and Ln i are the mass density and the length of an element, respectively; E and σ are the Modulus of elasticity and stress, respectively.
The first objective is the mass of the structures and the movement of the nodes, or connection points, throughout the truss construction, and is called nodal displacement, which is the secondary objective.We can ensure that the structure will deflect sufficiently under loads without exceeding the critical limits by optimizing the nodal displacement.
We used dynamic penalty functions to manage the constraints associated with compressive and tensile stresses [64].The penalty function penalizes the objective functions based on the extent of the constraint violation.The penalty function is formulated as follows: where, p i is the value of constraint violation concerning the bound p * i .C is a penalty function that calculates the total constraint violation, where Ci represents the extent of the breach for each constraint.The penalty function adjusts the objective functions based on the total constraint violation, with a penalty factor determined by constants ε 1 and ε 2 .The constants values, ε 1 and ε 2 are taken three based on comparable research findings from the literature [65].The primary objective is to minimize the maximum nodal displacement and total mass by optimizing design variables using Finite Element Analysis (FEA).Stress and cross-sectional area constraints must be followed during this optimization process to assure structural integrity and safety.

Truss structural problems
Five standard benchmark truss problems-the 10-bar, 25-bar, 60-bar, 72-bar, and 942-bar trusses-are used to evaluate the performance of the studied algorithms.Furthermore, the published results and those from previous studies are compared [56].The following sections provide details of these five common benchmark truss problems.5 shows the density and Yong's modulus of the truss material with loading conditions at different nodes.
The third benchmark truss is a Ring truss with 60 bars.[56,[66][67][68], as seen in Fig 5 .Table 6 illustrates how the truss is divided into 25 parts, each representing a symmetry.This truss is subjected to three load cases with size variable (cross-sections) values ranging from 0.5 in2 to 4.9 in2.Material properties are displayed in Table 6.
The 4 th benchmark, which is a 72-bar 3-D truss [56,[66][67][68], is presented in Fig 6 .As shown in Table 7, the truss is made up of 72 elemental cross-sections that are arranged into 16 segments.Two different load cases are considered as per Table 7.The maximum allowable stress is 25 ksi as a constraint to prevent structure failure.Material properties are as per Table 7 for the 3-D 72-bar truss structures.Here, 72-bar elements are grouped into 16 members as design variables.Size variables represent possible cross-sectional areas that can be assigned to truss structures ranging from 0.1 in 2 to 2.5 in 2 .
As seen in Fig 7, the fifth and last truss is a large tower truss with 942 bars [44,56,67] To preserve symmetry, the 942 members of the truss geometry are divided into 52 members, representing design variables as per Table 8.Two loading conditions are given: nodal conditions and lateral loading, which are applied at different truss portions.Size variables for this huge tower range from 1 in2 to 200 in2.Material properties are density and Young's modulus, as per Table 8.
This study rigorously tested six different algorithms, with each approach undergoing one hundred runs for each studied truss problem.These evaluations were conducted meticulously, with each test encompassing 50,000 functional evaluations, ensuring a comprehensive analysis of algorithmic performance across various scenarios and challenges inherent in truss design problems.

Empirical assessment
MOHO was applied to evaluate all considered truss problems and ascertain the effectiveness of the approximate Pareto-optimal solutions generated by the MO optimization algorithm.The ▪ The hyper-volume (HV) index measures the percentage of target space that members of the ND solution occupy set S. It provides information about the S set's convergence and diversity.A hypercube vi is created for every solution i in S by a collection of reference points.A higher HV value denotes an algorithm that performs better.A visual representation of HV for MO problems is shown in Fig 8.
▪ The gap between the actual Pareto-optimal front and the estimated Pareto-optimal front found during the search process is measured by the Generational Distance (GD).The actual and approximate Pareto-optimal fronts coincide when the GD value is zero.

IGD ¼
ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi In the given equation, |P| denotes the count of outcomes in the Pareto front, where di represents the Euclidean distance across the nearest solution from the reference front and the objective function vector of the ith solution in the acquired front.Conversely, |P'| signifies the quantity of the solutions for the reference plane.This metric is utilized to measure both front expansions and progression.Pictorial representation of GD and IGD matrices are shown in Fig 9 which assesses the proximity of solutions in the approximation set to the true Pareto front.
▪ Together, the metrics for extent (ET) and spacing (SP) create a new evaluation matrix called the spacing to extent (STE) ratio.This matrix allows one to examine the extent and spacing aspects simultaneously.A more efficient and non-dominated Pareto front usually has a lower STE value.
The Euclidean distance di separates the objective function vector of the i th solution from its closest neighbor.� d is the mean value of all di, where M is the number of objective functions.f max

Results-discussions and comparative study
The outcomes of the algorithms that are being examined are shown in below Tables 9-12, which are elucidated as follows: Table 9 shows the hypervolume (HV) measures, which indicate how the algorithms' ND sorting capability has evolved and changed over time.Greater HV values are associated with superior non-dominated fronts.Out of all the problems examined, MOHO had the highest HV values, followed by MOMFO and DEMO, which suggests that it performed better in reaching a reasonable convergence rate.▪ Variation in Hypervolume (HV): For all five considered truss constructions, MOHO shows significantly better average, maximum, and minimum HV values than MOMFO, MOALO, DEMO, NSGA-2, MOACS, and MOAS.Furthermore, MOHO's standard (std) value from the HV test is excellent for all the algorithms considered for all truss constructions, indicating better search consistency.
▪ Performance Comparison: Among all the average HV values of the algorithms considered, MOHO exhibits the highest values for truss structures, with respective averages of 53518.79,1894.07,4087.08,2317.36, and 75288799 for 10-bar, 25-bar, 60-bar, 72-bar and 942-bar truss respectively.Algorithms with higher average HV values are more effective at exploring the solution space and identifying Pareto-optimal solutions.▪ MO algorithm effectiveness: Comparing the HV values for each MO algorithm across the different truss structures allows us to identify which algorithms consistently perform well across various optimization problems.
▪ Standard deviation analysis: Std value by the MOHO for all considered trusses are less compared with other optimizers with Friedman values of 1.10, 1, 1, 1.03, and 1 for 10-bar, 25-bar, 60-bar, 72-bar and 942-bar truss respectively.The stability and consistency of ▪ Sensitivity of truss structures: The MOHO framework balances exploitation, prioritizing feasible solutions, and exploration, which involves seeking diverse solution options, and offering a well-rounded optimization approach.With these capabilities, MOHO holds promise in systematically navigating the vast solution space of larger structures, potentially approaching optimal designs.Rigorous comparisons using Friedman's rank test rank algorithms, with MOHO emerging with the top score of 1.03 among all algorithms with conflicting objectives, indicating superior Pareto front quality compared to alternatives.Analyzing the HV values in the table enables a comprehensive assessment of multi-objective optimisation algorithms' effectiveness, robustness, and sensitivity across a range of truss structures, facilitating informed decision-making in algorithm selection and application.
Figs 11-15 depict the Pareto fronts for truss problems across seven algorithms.They reveal the correlation between mass and maximum displacement, which provides visual insights into the trade-offs between these objectives.Each figure represents a scatter plot showing the   ▪ Algorithm Performance: The distribution and spread of points along the Pareto front provide insights into the performance of each optimization algorithm.Algorithms that generate solutions distributed across a wide range of Pareto front regions demonstrate better exploration capabilities and versatility in identifying diverse trade-off solutions.
This observation underscores MOHO's superior performance.Its Pareto fronts are characterized by smoothness and even distribution, in stark contrast to the fragmented and discontinuous fronts generated by other multi-objective algorithms.This distinction highlights MOHO's prowess in optimizing truss structures more efficiently and effectively across various design objectives.
The hypervolume generation process with the function evaluation of different truss constructions is seen in Figs 16-20.The hypervolume metric's evolution during the optimization process is depicted graphically in these figures, which also highlight the various algorithms' abilities to explore and converge towards the Pareto front.By scrutinizing these figures, it becomes apparent how effectively each algorithm performs concerning coverage and convergence toward optimal trade-off solutions.Also, the figures illustrating hypervolume through function evaluations for all MO optimization algorithms across various truss structures provide insights into the performance of each algorithm over time.By plotting the hypervolume   10, which is an essential metric for evaluating the differences between the Pareto optimal front and ND solutions in various truss configurations.A lower GD score indicates an exceptional, non-dominated front.For the 10-bar truss average GD value is 4.1908 with a standard deviation (std) of 0.5805, for the 25-bar it is 0.4071 with a std 0.1666, for the 60-bar it is 0.6826 with a std of 0.3439, for the 72-bar it is 1.4976 with a std of 1.0042 and for the 942-bar truss structure average GD is 1609.70 with a std of 230.60 for MOHO.The GD measure results reveal that MOALO, MOACS, MOMFO, and MOHO have excellent non-dominated fronts and perform extraordinarily well, which means they can generate solutions closer to the true Pareto front and well-distributed across the objective space.On the other hand, MOAS, DEMO, and NSGA-2 perform worse when evaluated using the GD metric.These comparative analyses help to understand the performance of all considered MO algorithms for all considered five truss structures in terms of their ability to explore the solution space effectively and produce high-quality Pareto fronts.informed decisions with the help of these figures, which are crucial tools for comprehending the dynamic behavior of optimization algorithms.
Spacing represents a good distribution of solutions with evenly spaced points, indicating that the MO algorithm has explored the objective space effectively and generated the most diverse solution set.Extent is a range of solutions along the objective space; a higher value means the solutions cover a wide range of objectives, indicating diversity.By simultaneously evaluating both spacing and extent, the Spacing-To-Extent (STE) metric offers essential information about the quality of non-dominated fronts.A lower STE score indicates a superior, ND front because it matches the spacing-to-extent ratio better.Table 12 presents the STE findings for the MO algorithms under consideration in the study.For 10-bar, 25-bar, and 72-bar trusses, the average STE value of MOHO is 0.0110, less than MOAS, MOACS, MOALO, and NSGA-2.For a 60-bar ring truss, an average STE value from MOHO is 0.0077 with the least std of 0.0089, the third best compared with other MO algorithms.Same for the large tower truss 942-bar, an average STE measure is 0.0088 with a std of 0.0017 with a Friedman value of 2.3.As can be seen from the STE findings, DEMO performs best, demonstrating its ability to create comprehensive, evenly spaced, non-dominated fronts.As the runner-up, MOMFO and MOHO exhibit STE results similar to DEMO's.According to Friedman's test at a 95% significance threshold, while other metrics may yield comparable findings, the overall Friedman test positions DEMO, MOMFO, and MOHO as the top three algorithms, reaffirming their effectiveness in generating high-quality non-dominated fronts with favourable spacing to extent ratios.
Figs 26-30 depict swarm plots representing the performance of optimization algorithms across all the truss structures, respectively.Swarm plots visually represent the distribution of objective function values or performance metrics obtained by different algorithms for each truss structure.These swarm plots provide information about how the solutions produced by each algorithm are distributed among the various truss structures.This analysis aids in   An overview of the comprehensive Friedman rank for all truss structures achieved by the techniques under consideration is shown in Table 13.The average Friedman's score for MOHO is 2.2850, which is the lowest and has the first rank compared with MOMFO, DEMO, MOALO, MOACS, MOAS, and NSGA-2.MOHO has an excellent convergence rate compared with the other prominent MO optimization algorithms.This dominance of MOHO is statistically significant, as indicated by Friedman's rank test at a 95% level, further underscoring its better performance than other algorithms assessed in the research.Overall, MOHO has the most excellent HV values, indicating that it explores well and comes up with various solutions.The results of MOHO, MOALO, MOACS, and MOMFO are near the optimal Pareto front, as indicated by their lowest GD values.In most instances, MOHO also had the lowest IGD, indicating a decent mix of convergence and diversity.MOHO was the most effective algorithm for these truss structure problems when ranked according to all three measures.Put more, MOHO is the best preference for these challenging engineering investigations because it finds a good distribution of near-ideal, well-balanced solutions.

Conclusion
This article presents the MO version of the novel hippopotamus optimization algorithm for solving five structural truss problems.Minimization of both objective function, structural mass, and maximum nodal displacement is subject to stress and area constraints.This algorithm's two decisive exploration phases and one exploitation phase generate excellent results for the truss optimization problems to examine its exploratory, exploitative, local optima evasion, and convergence properties.Through quantitative and qualitative analyses, comparing MOHO with six prominent algorithms based on four significantly recognized performance measures, we demonstrated its effectiveness in handling real-world truss structure More research into how it performs on higher-dimensional engineering design challenges is essential to evaluate MOHO's potential thoroughly.Further studies could examine how MOHO can be applied to multi-modal and multi-dimensional technological problems with conflicting objectives.Additionally, the research can be expanded to investigate methods for enhancing performance and carrying out evaluations compared to other well-known optimization techniques.Further advances in efficiency and scalability will enable the MOHO for truss structures to manage more prominent and intricate structural systems.Promising directions include investigating its adaptability to different structural types beyond trusses, integration with advanced analysis techniques, and resilience in managing unknown parameters.The applicability and relevance of MOHO are further expanded in various engineering areas by   hybridization with complementary optimization approaches and application to real-world engineering challenges.Enhancements in the future could involve extending MOHO to handle more extensive and more intricate truss structures containing increased design variables and constraints and additionally, adapting to dynamic environmental changes in design requirements using adaptive parameters.Overall, MOHO exhibits potential as a proficient and successful method for multi-objective optimization in truss-bar design issues, and its broader applicability to optimization contexts demands more research and analysis.

Fig 2 .
Fig 2. Flow chart of the MOHO.https://doi.org/10.1371/journal.pone.0308474.g002 Figs 3-7 show the investigated truss constructions and their loading circumstances, with geometrical dimensions used for computational evaluations.Tables 4-8 summarize the truss design considerations.The first structural problem is a 10-bar truss, depicted in Fig 3. Angelo et al. have extensively studied this 2-D truss [56, 65-68].Ten design variables were used to solve MO problems.As depicted in Fig 3, a vertical downside load of 100 kips is applied at nodes 2 and 4.

Fig 4
Fig 4 shows the 25-bar truss used in the second truss problem.Academic research regularly utilizes this 3-D truss[56,[66][67][68].It is symmetric about the x-z and y-z planes, and its twentyfive elements are arranged into eight groups.Table5shows the density and Yong's modulus of the truss material with loading conditions at different nodes.The third benchmark truss is a Ring truss with 60 bars.[56,[66][67][68], as seen in Fig 5.Table6illustrates how the truss is divided into 25 parts, each representing a symmetry.This truss is subjected to three load cases with size variable (cross-sections) values ranging from 0.5 in2 to 4.9 in2.Material properties are displayed in Table6.The 4 th benchmark, which is a 72-bar 3-D truss[56,[66][67][68], is presented in Fig 6.As shown in Table7, the truss is made up of 72 elemental cross-sections that are arranged into 16 segments.Two different load cases are considered as per Table7.The maximum allowable stress is 25 ksi as a constraint to prevent structure failure.Material properties are as per Table7for the 3-D 72-bar truss structures.Here, 72-bar elements are grouped into 16 members as design

▪
distribution of solutions along the Pareto front, where each point represents a unique solution generated by the optimization algorithm.▪ Truss Structures: The figures display Pareto fronts for different truss structures, such as 10-bar, 25-bar, 60-bar, 72-bar, and 942-bar, allowing for comparisons across various design complexities and sizes.▪ MO Optimization Algorithms: The plots include Pareto fronts generated by seven multiobjective optimization algorithms, such as MOAS, MOACS, DEMO, NSGA-2, MOALO, MOMFO, and MOHO.Each algorithm employs distinct optimization strategies and techniques to explore the solution space and identify Pareto-optimal solutions.Mass vs. Displacement: The scatter plots' axes represent the two objectives of interest: mass (structural weight) and displacement (structural deflection).The position of each point on the plot indicates the corresponding values of mass and displacement for a particular design solution.▪Trade-Off Analysis: The researcher can analyse the trade-offs between both objectives by examining the distribution of points along the Pareto front.Solutions located closer to the Pareto front represent superior trade-offs, where improvements in one objective (e.g., reducing mass) come at the expense of the other (e.g., increasing displacement).

Table 3 . Stages of the exploitation phase.
t ; t ¼ 1; 2; . . .; T max t is the current iteration, and T is the maximum iteration.